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In this talk I present the logic behind, and examine the reliability of, estimates of the 
critical end point (CEP) of QCD using the Taylor expansion method. 

> . 

1. The method and results 

The global structure of the expected phase diagram of two flavour QCD is the following 
[Q]. In the chiral limit of two flavour QCD, i.e., when m n = 0, one expects the chiral 
phase transition to occur at a critical point at finite temperature, T, and zero baryon 



O 



chemical potential, [lb- This critical point is expected to develop into a critical line in the 
plane of T and hb for m n = 0, and turn into a line of first order transitions at a tri critical 



point. When ^ 0, the structure changes dramatically. There is no critical line; rather 
there is a first order line ending in a critical end point (the CEP) in the Ising universality 
class at T*(m„-) and [i B (m n ). With varying the CEP traces out a "wing critical line" 
ending at the tricritical point. These considerations are based on symmetry arguments 
and universality, and hence are expected to be robust. 

Direct verification of this picture involves lattice QCD computations in the chiral limit 
and finite chemical potential. Both are technically unfeasible at this time, the former due 
to chiral slowing down (computer time requirements diverge as m 7T — > 0) and the latter 
due to the sign problem (Monte Carlo is impossible at finite hb)- The strategy of [H] is 
not to prove or disprove this picture directly, but, by assuming this picture, to estimate 
T*(m 7r ) and ^(m^). 

This is done by making a Taylor expansion of the pressure around /ib = — 

P{T,Hb) = (|) IogZ(T^ fl ) = P(T,0) + T,Xb(T)^-. (1) 

CP symmetry of the problem reduces every odd coefficient to zero, hence the sum starts 
with the term in n = 2. This induces a series for the quark number susceptibility (QNS), 
which is the second derivative of the pressure with respect to fis, and the associated 
radius of convergence — 
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We have used T to convert the series to one in the dimensionless variable Hb/T. Since we 
use this expansion at fixed T, there is complete equivalence with the dimensional form. 
A notational convention used here is x(T,Pb = 0) = Xb(T)- 

At the CEP this QNS is expected to diverge in the thermodynamic limit. Let p* denote 
the limit of p n as n — > oo, evvaluated in the thermodynamic limit. Then, somewhere on 
the circle of radius p*(T) in the complex plane, we are assured of finding a singular point. 
If this singular point is on the real axis, then we have found the critical end point. 
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Figure 1. p B as a function of the order, Figure 2. p 5 ' 7 as a function of the lattice 
n, for small (Lm n ~ 3, circles) and large size, Lm v at T = 0.95T C . 
(Lm n ~ 9, boxes) volumes at T — 0.95T C . 



One complication is that in a lattice computation one perforce uses finite spatial 
volumes — a cube with sides L. On finite volumes there is no true phase transition, 
so the limit of p n as n — > oo must be infinite. One expects the sequence p n (T,L) to 
approach some value p*(T,L) for restricted values of n, before climbing to infinity. As L 
increases, the range of n for which this stable plateau in p n (T, L) is seen should increase 
and the plateau value, p*(T,L), should converge to the thermodynamic limit p*(T). In 
Figure ^ we show that p n (T,L = 3.3/m n ) increases continually with n, whereas a stable 
value is seen for p n (T,L = 10/m n ). In Figure 121 we show that p 5,7 (T,L) cross over from 
small to large volume behaviour at Lm n > 5. The extrapolation from large volumes to 
the thermodynamic limit can be made using finite size scaling. In [|2] this is performed 
assuming that the CEP is in the Ising universality class. 

In [ |2] as T was lowered from a very high value, it was found that the radius of con- 
vergence decreased from a very large value to p* ~ T. Oscillations in the sign of the 
non-linear susceptibilities (NLS, \b) indicated that the singularity was at complex p*. 
On further lowering T it was found that all the series coefficients have positive sign (Fig- 
ure EJ), showing that the singularity moves to real p*, allowing the identification of T*(m,r) 
and p B (m n ). 

The technical part of the computation involves the choice of action. For the gauge fields, 
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Figure 3. The temperature dependence of 
the successive NLS. 



Figure 4. Variation of p 3 with m,,-; the bar 
is data from [|7j. 



we have taken the Wilson action and for the quarks we have two flavours of staggered 
quarks. Staggered quarks at finite chemical potential have many uncertainties at finite 
lattice spacing [ El E] . Since the method of Taylor expansions uses only simulations at 
[Mb — [El Ej, they may give results which agree with those obtained by reweighting [E] 
or analytic continuation of simulations at imaginary /xg [ 6J only in the continuum limit. 
At present there are no simulations with Wilson quarks. 

The lattice simulations were carried out at lattice spacings a = 1/4T. The scale was 
set by the finite temperature crossover at T c , giving m p /T c = 5.4. The quark mass was 
chosen such that m n /m p = 0.31 (compared to 0.18 in the real world). This agrees with 
the simulation parameters of [ El IE] but not with those of [ [Zj which differ by having 
m^/rrip = 0.7. 

An exploratory study of the variation of the CEP with m n /m p can be performed in 
partially quenched QCD, i.e., in QCD where the sea quark masses are allowed to be 
different from the valence quark masses. In such an exercise we found the radius of 
convergence to vary by large amounts. We show our results in Figure EJ It is interesting 
that the results of [|7j are close to the curve generated by these computations, indicating 
that the apparent difference between and [|7] could largely be a matter of change in 

TThr. 



2. Two technical points 

A major question in lattice simulations is that of statistics. In [ [2] the measurements 
were taken on N^ t e / {2T mt + 1) = 50-200 statistically independent configurations. Sta- 
tistical independence is measured by the integrated autocorrelation time, T mt . Since sta- 
tistical errors change as J (2r mt + 1)/N^[ t e , increase of statistics by a factor of 4 could 
easily be offset if r int = 1 instead of zero. Much more important is the statistics used in 
the estimator of fermion traces, which we shall discuss next. 
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2.1. Evaluating traces and their products 

The Taylor expansion coefficients, the NLS xS; m e Q s - (DEI) involve derivatives of the 
quark determinant, which become traces over quark loops. Traces of very large matrices 
are efficiently evaluated by the identity Tr A = Tr AR where R is a stochastic represen- 
tation of the identity matrix — 

i N v 

R=^rY,\r>)(r>\ N ^°I, (3) 
ZiV *> i=i 

the limit being obtained when |rA is one of a set of N v vectors of complex numbers, each 
component drawn independently from a distribution with vanishing mean. Typically one 
has to evaluate traces as well as products of traces — 

a j = 2 ^2^2 A ^ r C r ip, Pn = Y[a j , p^ = Y[TrA j = H 

a/3 t=l j=l j j 

where the factorisation of the average (denoted by a bar over the symbol) of p n occurs 
when the different sets of random vectors r J are independent. Here a and (3 are matrix 
indices. If the random vectors are drawn from a Gaussian distribution, then a 3 - is clearly 
a Gaussian distributed variable with mean a] = Tr Aj and variance <rj oc 1/N V . Hence 
the product of traces is a product of independent Gaussian random variables. 

A product of Gaussian random variables is not Gaussian. To prove this, consider 
8j = aj — aj, which is Gaussian distributed with zero mean, so = cr| and (Sj) = 3a j, 
so that [6j] = (8j) —3(Sj) = 0, as must be true for a Gaussian. In that case, however, the 
product d = n?=i $j nas the property that 

n n n n n 

(d 2 ) = n ® = n °i (d 4 ) = n = 3n n ^ = ( 3n - 3 ) n ^ ^ o ' ( 5 ) 

j=l 3=1 3=1 3=1 3=1 

where the factorisation of the expectation values follows from the independence of the 
6/s. Thus products of Gaussians are highly non-Gaussian. 

Since the central limit theorem applies, estimators of d are still Gaussian distributed. It 
would seem that the fourth cumulant varies as (3" — 3)/A^ n , and therefore the Gaussian is 
approached rapidly. In actuality the situation is worse, as one can check by reconstructing 
the distribution of p n = Ylj(fy J r Sj). We find that [p^J is dominated by the slowest falling 
term, which varies as n 2 /N^, rather than the most rapidly falling term (3/A^) n . It was 
found that N v ~ 30 are needed to control errors for n = 2 [ 3j. Consistent with this, 
N v ~ 500 are needed to control errors for n = 8 [Ej. Conversely, using N v = 100 for 
n = 6 [Ej is equivalent to working with ~ 10 for n — 2. 

2.2. Critical divergence 

Since we analyze the series expansion of eq. (0) to estimate the critical end point 
shouldn't one be able to use the series to see critical divergence? In fact, the lack of such 
a divergence has been used as an argument for the absence of the critical end point in [Cj. 
However, it has been known since the 1960's that the partial sum of the series expansion 
is generically perfectly smooth across the critical point as determined by the radius of 
convergence upto the same order. 
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Figure 5. Partial series for the ratio m/x- 
There is no hint of critical behaviour at 
y = l. 



Figure 6. A comparison of the CEP ob- 
tained in [ |2] with the freezeout curve in 
heavy-ion collisions [Ej. 



As an illustration I show results for the series expansion of the magnetic susceptibility, 
X, of the two dimensional Ising model on a square lattice. The expansion is made in 
y = exp(—2H/T) around the point y = 0, where H is the magnetic field and T the 
temperature. This series has a radius of convergence unity. The exact ratio m/x (where 
m is the magnetization) vanishes in the ordered phase. However, as shown in Figure 
the ratio is perfectly smooth across y = 0. 

One consequence of this realization is that the agreement of a series expansion of the 
pressure summed to low orders with any effective theory, for example, a resonance gas 
model, cannot be used to exclude the possibility of a critical end point. 

3. Summary 

There are four measures of reliability of lattice simulations of QCD — 

1. the lattice spacing, a, which must be taken to zero, but in current computations is 
1/4T [ EH El • The next generation of computations will reduce the lattice 
spacing by 33%. 

2. the quark mass, which has to be tuned by fixing the pion mass or the ratio m n /m p . 
The physical value of this ratio is 0.18. In [0EIE1 this was taken be 0.31, in [|7j it 
is 0.7 and in [ [H] it is the physical value. 

3. the spatial size of the lattice, L, which must be taken to infinity through the process 
known as finite size scaling in order to obtain reliable results for thermodynamics. 
It was found that Lm n < 5 gave rise to tremendous finite size distortions [|2J. Only 
[ 13 Hj use lattices larger than this. 

4. the statistical precision, N^ t e = y / A 7 s b ^ e /(2r in4 + 1), where N^ t a J t e is the number 
of configurations generated and r int is a measure of statistical independence called 
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the integrated autocorrelation time. When T- mt = 0, the naive statistics actually 
represents the true statistics, but even the small value r- mt = 1, causes the true 
statistics to be three times smaller than N^ t e . In [ E| A£™ e = 50-200 at all T, 
whereas no other simulation reports measurements of T int . 

In addition, any computation of the expectation of products of fermion traces contains 
another measure of reliability. This is the number of vectors, N v , used in evaluating the 
traces. Its crucial role is explained in Section 12.11 

The method of Taylor expansions can give reliable answers to questions about the QCD 
phase diagram once all these technical matters are under control. The critical end point 
can be determined in two steps — 

1. Generate the series coefficients for the QNS in eq. (J2J) and evaluate the radius of 
convergence of the series. 

2. When all the series coefficients have the same sign, the radius of convergence implies 
a singularity at real yU#, which is the location of the CEP. 

This works although the partially summed series for the pressure (eq. or the QNS (eq. 
12} are not expected to give the physical values of these quantities, as we show through 
examples in section 12.21 As a result, comparison of models, such as the resonance gas 
model, with the partial series sums is blind to the physics of the critical end point. 

The CEP of QCD can be obtained from the results presented in [|2], with Lm^ > 5, pion 
mass within 50% of the physical pion mass, and with good statistical control. The result 
is shown in Figure HI superposed on the freezeout curve obtained in Computations at 
smaller lattice spacings would be very useful in checking the stability of this result. 
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